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A countable set of asymptotic space - localized solutions is constructed by the complex germ method 
in the adiabatic approximation for the nonstationary Gross - Pitaevskii equation with nonlocal non- 
• linearity and a quadratic potential. The asymptotic parameter is 1/T, where T 3> 1 is the adiabatic 

evolution time. 

A generalization of the Berry phase of the linear Schrodinger equation is formulated for the Gross- 
Pitaevskii equation. For the solutions constructed, the Berry phases are found in explicit form. 



> 

Introduction 

in 

' A quantum system with slowly (adiabatically) varying parameters is characterized by the conservation 
of its quantum numbers throughout the time of adiabatic evolution when the system Hamiltonian has a 
nondegenerate spectrum. In other words, the vector of the system state remains an eigenvector at any time 
f*^ , and gains a phase factor during the adiabatic evolution (see, e. g., [1]). M. Berry has shown [2]) that the 
phase gained consists of a dynamic and a geometric (topological) part. Such a division is caused by different 
mechanisms of their origination. The dynamic phase is associated with the evolution of the system and the 
topological one with the geometry of its parameter space. 

Both phases result from the decomposition of the general phase of the leading term in the semiclassical 
' approximation of the solution of the Schrodinger equation in a small parameter 1/T, where T is the adiabatic 
^ , evolution time. The dynamic phase corresponds to the zeroth order term in the 1/T-expansion, and the 

• Berry phase belongs to the first order term. Inclusion of the topological phase in the general phase is 
. 1-H , necessary to determine the leading term of the semiclassical asymptotic series. The work [2] has initiated 

• comprehensive investigations of the geometric phases of linear quantum- mechanical equations. For details 
^ ', see the reviews [3-6]. 

■ - - ' It is of interest to study geometric phases for nonlinear equations describing various class of nonlinear 
phenomena [7] whose nonlinear properties are due to collective interactions. An important example of 
systems showing nonlinear properties is the Bose-Einstein condensate (BEC) [8,9]. 

In the BEC theory various models are used. In [10, 11] the BEC model is based on the two-level Hamil- 
tonian of a system of particles confined by a magnetic field. The system has a constant number of particles 
distributed between two subsystems, which are in different states. The Berry phase for this system is ob- 
tained under the condition that the system parameters responsible for the interaction of the subsystems vary 
slowly. 

Many of the BEC models are based on the Gross-Pitaevskii equation (CPE) [8,9, 12, 13], which is a non- 
stationary multidimensional cubic-nonlinear Schrodinger equation with the potential of an external field. 
The BEC states are described by localized solutions of the GPE. In studying these models, a number of 
serious mathematical problems arise. We mention two of them which are essential for our consideration. 

First, no method of exact integration for the GPE are available except for the nonstationary one- 
dimensional Schrodinger equation that is integrable by the Inverse Scattering Transform method [14] with 
no external field. In the presence of an external field, only an approximate solution can be found in terms 
of the soliton perturbation theory [15] with the assumption of a weak external field. 

Second, even for a two-dimensional space, the localized solutions of the nonlinear Schrodinger equation 
(NLSE) with focusing cubic nonlinearity and no external field are unstable and eventually collapse [16], 
which is not observed in experiments. On the other hand, models based on the GPE with local nonlinearity 
can be considered as simplified versions of models with the GPE having a nonlocal nonlinearity [18]. This is a 
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reason to consider in more detail the nonlocal operator J V{x — t)\'^dy that arises in the derivation 

— oo 

of the Gross-Pitaevkii equation. The interatomic potential V{x — y) is usually assumed to be short-range, 

and therefore the above nonlocal nonlinear operator can be replaced by a local operator /3|^'(x, t)p, where 
+00 

P = J V{y)dy. More detailed analysis of the nonlocal properties can be found, for example, in [18] where 

—00 

the wave function of the one-dimensional equation is expanded in a series in z = a; — y. The symmetry of the 

potential and its decrease at infinity result in an additional term oc -^-^ \ ip in the local Gross-Pitaevskii 
equation. The equations obtained are studied by numerical methods as the construction of an analytical 
solution fails even in the one-dimensional case. 

In this work we consider the one-dimensional nonlocal Cross-Pitaevskii equation and use an approach 
where the potential V{x,y) is expanded in the variables {x,y) to the second-order terms. Exact solutions 
are constructed for this equation in [19,20]. Although the interatomic quadratic potential does not decrease 
at infinity, the convergence of the integral is provided by a proper choice of the class of functions in which 
we seek solutions to the equation. 

The one-dimensional GPE describes BEG states if the longitudinal dimension of the condensate is much 
greater than its cross dimension [18]. Unlike the NLSE with local nonlinearity, the basic properties of the 
solutions of the nonlocal one-dimensional GPE conserve in the multidimensional case. 

The aim of this work is to obtain the Berry phase in explicit form for a one-dimensional GPE of the form 

ihdt + HUR{t),'f{t))^'f = 0, (0.1) 

n^{R{t),^{t)) = n{R{t)) + KV{R{t),^{t)), (0.2) 
n{R{t)) - ^^^^^^ , cr(Oa^^ , p{t){xp + px) 
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ViR{t),^{t))^- I dy ait)x^ + 2bit)xy + c{t)y' \^{y,t)f 

which involves a quadratic potential and an external field of the harmonic oscillator type. Here a(t), b{t), and 
c{t) are the potential parameters; >c is the nonlinearity parameter; and /i(t), cr(i), and p{t) are time-dependent 
system parameters. Therefore, the Hamiltonian depends on time via the set of parameters R(t)~{ii{t), crit), 
p{t), a{t), b{t), and c(t)). The quadratic potential in Eqs. HO.lfl . (|0.2|l models the magnetic traps that confine 
the condensate [17]. 

The problem of constructing solutions for Eq. H0.1|) with an external potential is solved with the use of 
the semielassical integration method developed in [19-22], where Eq. (|0.1|) is called the "Hartree type equa- 
tion" . Below we call the nonlinear operator in Eq. (|0.1|l the nonlinear Hamiltonian. The semielassical 
integration method allows one to construct localized solutions, approximate for the potential of general form 
and exact for the quadratic potential. We use these solutions to find and study geometric Berry phases. 

Let us give necessary facts from the theory of geometric phases of linear quantum mechanical sys- 
tems to find Berry phase in the nonlinear case. Following [2] (see also [3,5]), we consider a Hamiltonian 
7i(i) = 7i(i?(t)), which depends on time t via a set of slowly varying T-periodic functions R{t). Denote by 
E„(R(t)){Xi R{t)) the eigenfunctions of the instantaneous Hamiltonian 7i(i?(t)): 

H(i?(i))*£„(fl(t))(x,i?(i)) ==^„(i?(t))^'£;„(^(t))(a;,i?(t)). (0.3) 

Assume the spectrum of the Hamiltonian 7i(i?(t)) to be nondegenerate at any fixed time and set the Gauchy 
problem 

{-ihdt+n{R{t))}-^ = 0, (0.4) 

*lt=t„ = *B„(i?^(to))(a^'^(io)). (0.5) 

According to the adiabatic theorem [1],^ the quantum numbers of the system will conserve throughout the 
adiabatic evolution time T if the parameters R{t) depend on time adiabatically. In other words, a solution 

*It is generally agreed that a system evolves adiabatically if 

■ T 

max Rj — < 1, (0.6) 

i— l,n Ri 

where Ri are the parameters of the Hamiltonian (see [23]). 



of the Cauchy problem (|0.4|l . (|0.5|) differs from the initial state by the phase factor 



'fix, T) = cxp[i0„(T)]*£„(fl(T)) (x, R{T)). (0.7) 
Following [2], represent the phase 4>n{T) as 

,^„(T) = <5„(r) + 7„(r), (0.8) 

where 6n{T) is the dynamic phase given by 

T 

Sn{T) ^ J E„{R{t))dt. (0.9) 



The phase ^n{T) is called an adiabatic Berry phase, and for a linear Schrodinger equation it is determined 
by the expression 

7„(t) = i{^E,ARit))ix,RitmE„iRit))ix,Rm (0.10) 

or 

T 

IniT) = I J {^'E^(R(t)){x,RitmEURm){x,Rmdt ^ j^A'jdW. (0.11) 
c 

Here 

^"-(*j||>- (0.12) 

and C is a closed contour in the parameter space. The functions A" act as components of an "induced gauge 
field" since under the gauge transformations 

^n^CXp{i^„{R))'ifn (0.13) 

the quantities A" are transformed as 

^ + H- ('^■^^^ 

Expression H0.11() does not depend on the transformation H0.14|l as the contour C is closed. The dynamic 
phase 5n{T) characterizes the mean value of the system energy, and the geometric phase "(n{T) does not 
depend on the system dynamics. Berry's phase depends on the geometry of the system parameter space and 
on the type of the contour C. 

Assume the adiabatic theorem to be hold for nonlinear equations and define Berry's phase by relations 
(I0.8|l and H0.9|) . Formula (|0.11|) . being equivalent to (|0.8|) in the linear case, needs additional substantiation 
for nonlinear equations. 

In this work we use the approach developed in [24. 25] for linear equations to find Berry's phase. The 
approach is based on the exact (or approximate) solution of the Cauchy problem H0.4|) . H0.5|) which is 
expanded in an adiabatic parameter. For the linear equation (|0.1() the exact solution of the Cauchy problem 
is constructed using the method proposed in [19]. Here, the initial Hamiltonian in the initial condition of 
the Cauchy problem (|0.5|l is replaced by (|0.2|l and the solution is determined by two auxiliary systems of 
ordinary differential equations: the Hamilton-Ehrenfest system and the system in variations. Solutions of 
these systems are unknown when the coefficients are arbitrary functions of time. When the coefficients depend 
on time adiabatically, we can seek the solution in the form of an expansion in the adiabaticity parameter 
which is taken as l/T, where T is a "long" characteristic time, for example, the adiabatic evolution period 
of the system. 

If the Hamilton-Ehrenfest system and the system in variations are solved accurate to O ( 1 /T) , the equation 
considered is solved with the same accuracy, and thus we obtain a solution of equation 1)0. l|l in the adiabatic 
approximation. For such a solution. Berry's phase can be found in explicit form. 

1 The Hamilton-Ehrenfest system 

For a linear operator A we define its mean value in a state ^{t) as 

im) ^ ^^(vl;(i)|l|vl;(t)) = A^{t, h). (1.1) 



On the solutions ^'(i) of equation HU.4|I we have 



d{A{t)) _ /dA{t)\ I 



dt 



{^) + l{[HUt,^{t)),Am, (1.2) 



where [A, B] — AB — BA is the commutator of linear operators A and B. 

Similar to the linear case, we call H1.2(l the Ehrenfest equation. From this equation with A = 1 it follows, 
in particular, that the norm of a solution of equation l|0.4|l is conserved, that is, 

\\^{x,t)r = \\'f{x,o)r = m'- 

So it is convenient to use the nonlinear parameter iir = >f||^'jp instead of >f. 
Denote by 



+00 

= ^ / **(y,0{(Ap,)'(Ay)n*(2/,t)d2/, k,l = o;^, 

—00 

the moments of the (k + l) - tli order centered relative to x^{t, ft) and p*(i, ft). Here Apy = —ihdy —pqi(t, ft), 
and {(Apj^)'(Aj/)'"'} is a Weyl-ordered operator with the symbol {ApyY (Ay)'' . In particular, 

CTxxit, ft) = a''^'^\t, ft), 

are the variances of the coordinates and momenta and of the correlation function of the coordinates and 
momenta, respectively. 

Consider the first-order operators p and x and the centered second-order operators (Ax)^, (Ap)^, 
{AxAp + ApAa;)/2, where Ap = p — p*(i, ft). Ax = x — a;*(i, ft), p*(i, ft) = (p), and a;*(t, ft) = (.t). 
The Ehrenfest system for the mean values of these operators reads 

p = -ao{t)x - p{t)p, 
X = ^i{t)p + p{t)x, 

Oxx = 2ii{t)axp + 2p{t)axx, (1-3) 
O'xp = p{t)<Jpp — (y{t)axx, 
^pp — ~2p{t)app — 2a{t)axp- 

Here 

cro(i) = <j{t) + k{a{t) + b{t)), d{t) = a{t) + xa{t). 
For system (|1.3|l . let us set a Cauchy problem with initial conditions 



p\t=s = x\t=s = a;o, 

(7pp\t—s — ^Oppi ^xp\t—s — ^Oxp-} ^xx\t—s — ^Ox 



(1.4) 



We call (IO|) the second order Hamilton- Ehrenfest system (HES) related to equation 

Consider the HES H1.3|l as a dynamical system which is not related to equation (|0.4(l . Apparently, not 
all solutions of the HES can be obtained as mean values of the corresponding operators on the solutions of 
equation HU.4f) . For example, the mean values must satisfy the Schrodinger uncertainty relation 

2 ft2 . , 

(TppCTxx - cr^p> — (1.5) 

for the second-order moments (for the higher-order relations see [26]). It can readily be seen that the HES 
admits the trivial solution p = 0, x = 0, a^'^''^ = 0, k + I = 2. The left-hand side of (|1.5|) is the integral of 
motion of the HES (|1.3|) (see [27]). Hence, it suffices that the uncertainty relation be fulfilled at the initial 
time. The uncertainty relations will be fulfilled automatically if the initial conditions for (|1.3|) are taken as 

P|t=s = PO = PV'(^)' X\t=s ^ Xo ^ x^{h), 

<^pp f=s = a,/, (ft), (^xplt^s = a:,, '(ft), (7xx\t=s = (ft), 



where 'ip{x,h) is the mitial condition for equation 1U.4|I : 



^{x,h,t)l^^ = i,{x,h). (1.7) 

Denote a trajectory in an extended phase space by g = g{t, €.) G M^, where 

fl(t,£) = (p(t,G:),x(i,£),(Tpp(t,G:),(Tp,(t,G:),a,,(i,e:))\ (i.s) 

€={CuC2,C3,C4,C5y 

is the general solution of the Hamilton-Ehrenfest system (|1.3() and g is the operator column 

0= {p,x,{Ap)^,^{ApAx- AxAp),{Ax)^y . (1.9) 

Here C;, Z = 1, 5 arc arbitrary constants which can be expressed in terms of the initial conditions H1.4|l . The 
matrix B'^ is transposed to the matrix B. The system (11.311 can be rewritten in the form 



where 







0-%, 


3\t=s 


= 00, 








-Pit) 



















Pit) 





















-2p{t) 


~2d{t) 















Pit) 





~a{t) 




V 











Mt) 


2p{t) 


) 



2 The associated linear Schrodinger equation 

Let us seek a solution to equation (|0.1|) in the form of the ansatz 



^{x,t,h) = ip(^^,t,\/Pj exp ^(^^(t, £) + P(t,e:)Aa 



(2.1) 



Here the function ip{£,, t, ^/h) belongs to the Schwartz space § in the variable ^ = Ax/^/h and depends regu- 
larly on V^; Ax ^ X- X{t, £). The real functions S{t, €) and Z{t, €.) = {P{t, €.), X{t, £)) that characterize 
the solution are to be determined. 

Expand the operators in equation H0.1|l in a Tailor series in Ax = x — a;^(i, ?i), Ay = y — x$(i, h), and 
Ap = p ~ pq, {t, h) . Then equation (|().l(l takes the form 



{^ihdt + Mt, *) + {^zit, Az) + -{Az,fi,,{t, *)Az)}* = 0, 



io(i,*) 



p{t)pl{t,h) a{t)xl{t,h) 



p{t)x^{t, h)p^{t, h) 



-^caf^\t, h) + ^{a + 26+ c)4(t, h), 



^zit,^) = 
^zzit,-^) 



ii{t)p^{t, h) + p{t)x^[t, h) 

a{t)xis/ {t, h) + pit)pq! {t, h) + x{a + b)xq, (t, K) 

' Pit) Pit) \ 
, Pit) m i ■ 



(2.2) 



Let us associate the nonlinear equation H2.2() with the linear equation that is obtained from (|2.2() by formal 
substitution of the solution of the HES p.3|) instead of the corresponding mean values of the coordinate and 



momenta operators and second-order centered moments. The resulting linear equation is 

{^thdt +S}{t, £) + {^,it, £), Az) + i(Az, 9j,,{t, g:)Az)}$ = 0, (2.3) 

+ |ca(t),,(t, £, + |(a + 2b + c)X^{t, £), 



ii{t)P{tx) + p{t)x{tx) 

' f4t) Pit) \ 
. Pit) ^t) ] ■ 



We call Eq. (|2.3() the associated linear Schrddinger equation. 
By direct check we see that the function 

$o(x, t, £) -I 0, t, €) = cxp 1 1 (^(t, £) + £) Ax- + ^^^x^) } (2.4) 

is a solution of equation H2.3|l . Here 

t 

5(t,e:) = J {Pit,€)Xit,€)-MtX))dt; (2.5) 



and -B(t) and C(t) denote, respectively, the momentum and the coordinate part of the solution 

am ^ ( ^ ) , (2^6) 

of the system in variations 

a = J532z(t)a, a \t=s= oq, (2.7) 

related to equation (|2.3() . 

The normalizing condition H^-H^ = 1 yields Nn = (7rft)-i/4(| C(0) D'^/^^ 
Let us introduce the notation 

a(t) = NaiC{t)Ap - Bit) Ax). 

If C(i) and Bit) arc solutions of equations H2.7|l . then the operator a(t) commutes with the operator of the 
associated equation H2.3|l . Therefore, the function 



nix,t,€) -^{a+it)j $o(.T,t,e:), 71 = 0,00 



will be a solution of the Schrodinger equation (|2.3|) . Commutating the operators a~^it) with the operator 
of multiplication by the function $o(a;,i, £) =| 0, i, £), we obtain the following representation for the Fock 
basis of solutions of the linear equation H2.3|l : 



a 2 



= 4^ (k) "$o(x, t, €) (^-^) "(7^)" exp (-in Arg C(t))i?„ (^Ax^ , 



where i?n(?) are the Hermitian polynomials and Qit) ~ Bit)C ^(t). Finding Na = (l/V^hj cxp [— lArgC(O)] 
from the condition [a(i),a+(t)] = 1, we have 

$„(x,t,e:) =1 n,t,G:) = ^[a+(t)]" I o,t,e:) = 
1 



^ I 0,t)(^^j expH^(ArgC(0-ArgC(0))]i/„(O, 
I^A.. (2.8) 



Using the properties of Hcrmitian polynomials, we obtain the mean values of the momentum and coordinate 
operators and the corresponding variances: 

2;$„ ~ 0, = 0, 

(2,0) ^ g(t)(2n+l) 

(I'l) M ^ p(t)(27i + l) 

(0^2),, _^ _ , (2n+l) 

'2ImQ(t)' 



''(t, ft) = Cra;a:(t, ft) 



The functions ^n{x, t, £) are solutions of equation for properly chosen £, such that the solutions of 
the Hamilton- Ehrenfest system solutions H1.3|l coincide with equations (|2.9|l . Denoting this set of parameters 
by £„, we obtain 

*„(a;,i) = *n(a;,t,£„). (2.10) 
The subscript n in £„ implies that every function '^nix,t) has its own set of parameters £„. 

3 Eigenfunctions of the instantaneous nonlinear Hamiltonian 

To construct the Berry phase using formulas (|0.7|l and (|0.8|) , we have to solve a spectral problem for the 
instantaneous Hamiltonian (|0.2() in the class of functions H2.1() 

H^{R,1pn)lpn = En^pn■ (3.1) 

To solve this problem, consider the nonstationary Schrodingcr equation 

{~ihdt+H^{R,'i'{t))}^ ^0. (3.2) 
Solutions of equation (|3.2() of the form 

*(x,t) = cxp{ - ^En{R)t}Mx,R) (3.3) 

provide a solution of the spectral problem (|3.1|l where ij^nix, R) and En{R) arc the instantaneous eigenfunc- 
tions and eigenvalues of the Hamiltonian Ti.^(R,ipniR)), respectively. 

The Hamilton-Ehrenfest system (|1.3f) for the first-order moments related to equation (|3.2fl takes the form 

(p = -a,.-pp, ^^^^ 
I X — ^j.p + px, 



and for the second-order moments we have 



^xp ~ l-^'^pp ^(^xxj (3-5) 

^pp '^P^pp ^(TfTjjp, 



where 

CTo = o" + if(a + &), a = (7 + xa. 

Let us denote 



f] = aop - p"^, n = \/ap~ p^. (3.6) 

The spectral problem is associated only with the time-localized solutions of the system (|3.4|l , (|3.5|1 which are 
stable in the linear approximation. The localization condition holds when 

= gqp — > 0, = dp — > 0. 

In this case the general solution of the system (13.41) is given by 

x{t) = Ci sin nt + C2 cos nt, 

p{t) = - inci ~ PC2) cosCit- -(nc2 + pCi) sinnt. (3.7) 



Accordingly, for the system 1)3. 5|) we have 

ly-xx {t) = C3 sin 2^t + Ca cos 217t + C5 , 

cr^p(t) = - - pC4 ) COS 2m - - (vLCi + pCs) sin 2^ - -C5, (3.8) 



opp{t) = \{{p^ - r2^)C3 + 2pnCi) sin 2m + ^((p^ - ii^)Ci - 2pnc3) COS 2m + -C5. 

M M M 

Here C;, Z = 1, 5 are arbitrary constants. 

FoUowing equation l|1.8|l . denote the general solution of the Hamilton-Ehrenfest system 1)3. 4|l . 1)3. 5|l by 
g = g{t, £) G K^. To solve the spectral problem, we have need for the stationary solution of the Hamilton- 
Ehrenfest system ()3.4|l . (|3.5(l g = g{t,<ts) G that is obtained with the parameters taken as follows: 

= (0,0,0,0,C5)T. (3.9) 
The system in variations H2.7|l for the corresponding associated linear equation (|2.3|) becomes 

a{t) = ( a{t). (3.10) 



p p 

Let us set a Floquet problem [28] for the system in variations ()3.1()(l : 

a{t + T) = e'^^'^a{t). (3.11) 

The quasipcriodicity condition (|3.11() for solutions of the system in variations (|3.10() is sufRcient for the 
solutions of the Hamilton-Ehrenfest system in the linear approximation to be stable. A solution of the 
Floquet problem l|3.1()(l and H3.11)l . with the normalization condition 

{a{t), a*{t)} = 2i, {ai, 02} = (ai, J*a2), 

can be written in the form 

«W = ^f-^ + ^M. (3.12) 



Let us seek solutions to equation (|3.2)l in the form (|2.1() . Then the solution 1)2.4(1 of the correspondent 
associated equation 1)2.3)1 is given by 



Accordingly, for the Fock basis 1)2.8)) we obtain 

£) = exp { - znm} (^— j HnU—Axj^o{x,t,q. (3.14) 
Taking the parameters of the solution of the Hamilton-Ehrenfest system in the form 

£„ = (Ci,C2,C3,C4,C5)' = (0,0,0,0,^ ^^^^^^^^^ , (3.15) 
we obtain that that solutions of the associated equation 1)2.3)1 will be solutions of the original equation ()3.2)l : 



2h p J \y hp 

Hence, the eigenfunctions of the operator ()0.2)) have the form 

1 \i/4 /Six 1/4 r i p ^ I Q. „ f I n 



and the corresponding eigenfunctions are given by 



1 \ / kcp 



4 Adiabatic approximation 



Assume the evolution of a system is adiabatic. To this case (see (|(J.6|l l the parameters of the Hamiltonian 
slowly vary in time, and we can introduce, along with the "fast" time t appearing in the time derivative, a 
"slow" time s on which the parameters of the Hamiltonian {R{t)=R{s)) depend. 

Let the "fast" and the "slow" times be related as 



t 



where T is the evolution period of the system. 

For equation HO.lf) we consider the Cauchy problem 



{-z;i5t+7^^(i?(s),*(t))}* = 0, 
*(a:,t)|^^Q =^„(x,i?(0)). 



(4.1) 



(4.2) 
(4.3) 



where '0„(x,i?(O)) are the cigcnfunctions of the instantaneous Hamiltonian 'H^{R{0),'^{t)). 

As noted above, to solve equation (|4.2|l in adiabatic approximation wc have to solve the Hamilton- 
Ehrenfest system and the system in variations accurate to the first order of 1/T. 

The Hamilton-Ehrcnfcst system for the first-order moments can be written as 



-p' = -cro(s)a; - p{s)p, 
^x' = n{s)p + p{s)x, 



(4.4) 



and for the second-order moments wc have 

1 , 



2p{s)(Jxp + 2p(s)cr 



^'xp = p{s)crpp - (T{s)axx, 



T 

rji^pp 



(4.5) 



-2p{s)app - 2a{s)a^ 



where a' ~ da/ds. 

Let us seek a solution to the Hamilton-Ehrenfest system H4.4|l . (|4.5|l in the form 

xit) = x("'(^) + ^^^'H^l Pit) = P^'H^) + ^p^'Hs), m = + ^s(i)(.), 

where 



(4.6) 




Substituting H4.6|l in the Hamilton-Ehrcnfcst system H4.4|l . H4.5|l and equating the term of the same powers 
in l/T, we obtain 



0(s) 

V m J 



f 


1 








p(s) 








p{s) 




+ Ci 




a{s) 






V 




) 





( \ 

2p(s) VO(s)/ 
1 (p{s)^' 



(4.7) 
(4.8) 



The function a)^x (s) is determined by the condition of existence of the following approximation and has the 
form 

(P{s)\' , ^ (p{s)- 



4'J(s) = Cl 



f7(s) 



-21]3(s) V^i(s). 

A trajectory in the extended phase space g = g(t, £) G has the form 

0(i,e:) = (o,o,app(t,e:),ap,,(i,£),a^^(t,e:))\ (Cx,c^)\ 



(4.9) 



(4.10) 



Applying the change of variables H4.1|l to the system in variations H2.7|) we have 

1 



T 



a'{t) = Jf},,{t)a{t). 



Let us seek a semiclassical asymptotic solution to the system (|4.11ll as 



ait) 



„j(T<E.(6)+0(.s)) 



fit), 



(4.11) 



(4.12) 



/(t)=/(°)(,s) + l/W(,s). 
Substituting (|4.12|) in H4.11|l and equating the terms of the same order in 1/T we obtain 



(4.13) 



f^'Hs) = 



/'p{s)\' 



(P[s)V 



2n{s) \fi{s 

-p{s) + in{s) 
fi{s) 



Let us expand the vector /'•^^(s) in the basis vectors /^°-'(s) and (s) to obtain 

/«(.)= a(s)/(0)(s)+/3(s)/(«)*(s). 

Then 

^{s) ( p{s) — i^l{s)\' 



f3{s) 



V u(s) ) 



41]2(s) V p[s) 



(4.14) 



(4.15) 



(4.16) 



The function a{s) is to be determined from the condition of existence of the next-order approximation. We 
do not give the explicit form of this approximation, as it does not contribute to the terms of the order of 
l/T. 

The matrix Q{t) = B{t)C^^ (t) accurate to the first order in l/T reads 



The function 



$r(.,i,e:) = (/Xf^y/%,p 



nh \ I c{t) I 



^(ArgC(O)-ArgC(i) 



X cxp 



— ( - / xcis)a^^{t,€,h)dt + Q{t)Ax^ 



(4.18) 



is a solution of the associated equation H2.3|l accurate to the first order in l/T. Here (Jxx{t, £, fi) is determined 
from expressions (|4.8|l and (|4.9|l . and Q{t) is given by formula H4.17|l . For the Fock basis we obtain the 
following representation: 



(a:, *,£) = ^<i><,'')^exp[-zn(ArgC(i)-ArgC(0))li7„(O, 



'Im Q(o)(s) 



Ax. 



(4.19) 



In view of (|4.17(l . we define £„ as 



(4.20) 



Accordingly, a solution to equation H4.2(l is determined as 

^'^Hx,t) = ^^:'>{xM. (4.21) 
Note that if H4.3|l is taken as an initial condition for finding a solution to equation (|4.2|) . then 

^{x,t) = ^^^\x,t) + 0(^), (4.22) 



where 



exp 



-T I En{T)dT + i-fnis) flpnix, R{s)). 



The function (|4.23|) is a solution of equation (|4.2|l in the adiabatic approximation. The quantities 



E„{s) = hin 



>cc{s) 
2Im 



1\ f >cc{s)^{s) 
2Vl{s) 



are the eigenvalues of the instantaneous Hamiltonian Ti.:^(R(s) , ipniRis)) , and 7n(s) has the form 



ln{s)^-{n+-) / 0'(r)- 



kc{T) Im Q(i)(t) 



2 (ImQ(o)(T))2 



dT=- 
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"+2 



>irc(T) /i(T) 



2 f]2(^) 



(4.23) 



(4.24) 
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m:c{t) pl(t) ' 


M(r) 




2 f72(T)J 


20(r) 





(l)'{T)dT = 

)'dT, j7(t) = yRTyr^j^MM^y"^. (4.25) 



Considering the evolution during a period, we obtain 



*(")(x,T) 



exp 



T / £;„(s)(is + i7„(r) k(,°na;,o). 



(4.26) 



Using (|0.8|) and H0.9|) . we determine the dynamic phase 



(4.27) 



and the Berry phase 



7„(T)= (n+i 



1 - 



2 f]2(s) 



1 



20 



P 



— dp - -^dA* = <P Aldfi + A" dp. 



(4.28) 



c c 
The components of the "potential" A" in the parameter space are determined by the following relations: 



An 



An _ 



1 



1 



XC fl 

T?72 



LP. 



"+2 



XC n 

^ Tfi2 
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217' 



(4.29) 



The Berry phase can be presented in the form of an integral over a surface 11 in the space of parameters, 
supported by the contour C: 

7n(r)=(n+i) JJ L^{a{dpAdfi)+p{dfiAda)+fi{daAdp)}- 



xc 3/i 

2" 4175 
(7 = (T + xa 



{(T(dp A d/i) + p(d/i A da) + p.{dd A dp) } , 



(7 = (T + >ir(a + c), 11 = \/[o' + iirajp — p2 



(4.30) 



Expression (|4.23() confirms the statement of the adiabatic theorem for a nonlinear equation with the Hamil- 
tonian (|0.2|l in the class of trajectory-concentrated functions. 

The difference between the Berry phases for the nonlinear equation H0.1|l and for the linear Schrodinger 
equation [x = 0) is the variation of the frequency 17 and in the appearing additional summand 



XC p 1 
T 172 217 



dp — —dii] . 
U J 



(4.31) 



Note that the calculation of the Berry phase from the states H3.17|l by formula H0.11|l does not give in the 
term H4.31fl . Accordingly, expression H4.29|> for the "potentials" differs from that obtained by formula ()U.12|I 



by terms proportional to >cc. This is so because the states (|3.17|l do not satisfy the superposition principle. 
In the limiting case as iir ^ 0, the Berry phase coincides with the one obtained earlier in [29] (see also [24]). 

The classical analog to a Berry phase is known as a Hannay angle (see, e. g., [30]), a component in the 
dynamic part of the "angle" variable that arises in an adiabatically evolable integrable dynamical system 
described in terms of the "action ~ angle" variables. Hannay's angle has a nature similar to that of Berry's 
phase, and the relationship between Berry's phase 7 of a quantum system and Hannay's angle O of the 
corresponding classical system is given by the formula 

where / is the quantized action and n arc quantum numbers. The differentiation with respect to n is 
performed as if it would be a continuous parameter. In a nonlinear case, it is natural to associate the 
Hamiltonian with an analog of Hannay's angle by the formula H4.32|l 



c 

In the limit iir ^ 0, we obtain the well- known Hannay angle for the generalized harmonic oscillator [30] 



4.1 Conclusions 

In this paper we have constructed a solution for the one- dimensional nonstationary Gross-Pitacvskii 
equation H0.1|) and have found explicit expressions for the adiabatic Berry phase. The eigenfunctions of the 
instantaneous nonlinear Hamiltonian H0.2|l have been constructed by a scmiclassical method based on the 
Maslov complex germ theory [34,35]. This method, approximate in general, gives exact solutions for the 
spectral problem in the case under consideration. 

A classical analog of the Berry phase is the Hannay angle [30]. We have defined the Hannay angle in 
terms of quantum mechanics, since the nonlinear problem requires a special study of the "classical equations" 
corresponding to the nonlinear "quantum" Gross-Pitaevskii equation. In our consideration the role of these 
classical equations is played by the Hamilton-Ehrenfest system (II. 2f) . which has no Hamiltonian form relative 
to the standard Poisson bracket. 

The geometric potentials obtained H4.29|l . which determine the Berry phase, are Abelian. In the linear 
quantum mechanical case {x = in equation H0.1|l ) they can be treated as effective potentials of an electro- 
magnetic field. Such a situation takes place, e.g., in the theory of the molecular Aharonov-Bohm effect in the 
Born-Oppenheimer adiabatic approximation. The corrections arising due to this effect are experimentally 
established in chemical physics [6]. These potentials were also considered in [31]. 

In conclusion, we note that the non- Abelian Berry phase is of interest in quantum computations [32,33]. 
Such a phase appears for an instantaneous Hamiltonian having a degenerate spectrum, which arises in a 
multidimensional case. The approach developed in this paper admits multidimensional generalization and 
the results obtained can be used to study the relevant problems. 
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